SASE FEL with energy-chirped electron beam and its appUcation 
for generation of attosecond pulses |j 

. . . E.L. Saldin, E.A. Schneidmiller, and M.V. Yurkov 

O I Deutsches Elektronen- Synchrotron (DESY), Hamburg, Germany 

^ ■ Abstract 

^ , Influence of a linear energy cliirp in tlie electron beam on a SASE FEL operation is studied 

^D ' analytically and numerically using 1-D model. Explicit expressions for Green's functions and for 

cni : 

output power of a SASE FEL are obtained for high-gain linear regime in the limits of small and 
O , large energy chirp parameter. Saturation length and power versus energy chirp parameter are 
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calculated numerically. It is shown that the effect of linear energy chirp on FEL gain is equivalent 
to the linear undulator tapering (or linear energy variation along the undulator). A consequence of 
this fact is a possibility to perfectly compensate FEL gain degradation, caused by the energy chirp, 
(— I \ by means of the undulator tapering independently of the value of the energy chirp parameter. An 
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application of this effect for generation of attosecond pulses from a hard X-ray FEL is proposed. 
Strong energy modulation within a short slice of an electron bunch is produced by few-cycle optical 
laser pulse in a short undulator, placed in front of the main undulator. Gain degradation within 



.^ \ this slice is compensated by an appropriate undulator taper while the rest of the bunch suffers from 

\Q • this taper and does not lase. Three-dimensional simulations predict that short (200 attoseconds) 



high-power (up to 100 GW) pulses can be produced in Angstroem wavelength range with a high 
degree of contrast. A possibility to reduce pulse duration to sub-100 attosecond scale is discussed. 
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I. INTRODUCTION 

Start-to-end simulations [l| of the TESLA Test Facility Free Electron Laser (TTF FEL), 
Phase 1 j^ , have shown a presence of a strong energy chirp (energy-time correlation) within 
a short high-current leading peak in electron density distribution that has driven Self- 
Amplified Spontaneous Emission (SASE) FEL process. The energy chirp was accumulated 
due to the longitudinal space charge after compression. According to the simulations (that 
reproduced well the measured FEL properties), the energy chirp had a dramatical impact on 
SASE FEL saturation length and output characteristics. A similar effect takes place during 
the operation of VUV FEL at DESY in a "femtosecond mode" [3|, y, |5[. Such a mode of 
operation might also be possible in future X-ray SASE FELs. 

There also exists a concept of frequency-chirped SASE FELs (frequency chirp of SASE 
FEL radiation is correlated with energy chirp in the electron beam due to the FEL resonance 
condition) aiming at the shortening of radiation pulse with the help of a monochromator 
Energy chirp can also be used to tune the output frequency of an FEL with coherent 
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prebunching as it was demonstrated in the experiment at the DUV FEL facility [7|. Thus, 
a theoretical understanding of the energy chirp effect on the FEL performance is of crucial 
importance. 

Analytical studies on this subject were performed in |8] in the framework of one- 
dimensional approximation. The general form of a time-domain Green's function as an 
inverse Laplace transform was derived in [»]. It was then reduced to the explicit expression 
in the limit of small energy chirp parameter up to the first order, resulting in phase correc- 
tion (and ignoring the gain correction). This explicit solution for the Green's function was 
used to analyze statistical properties of a chirped SASE FEL in this limit. A second order 
correction to the FEL gain was presented in jq| but this result is incorrect. 

In this paper we study the impact of energy chirp on SASE FEL performance. We also 
find that FEL gain degradation can be perfectly compensated by undulator tapering. We 
discuss an application of the compensation effect for generation of attosecond pulses from 
X-ray FELs such as European XFEL J9| and Linac Coherent Light Source [lO|. 



II. GREEN'S FUNCTION 

Let us consider a planar undulator with the magnetic field 

H.,{z) = H^cos{27cz/X^) , 

where A^ is undulator period, and H„ is peak magnetic field. Electric field of the amplified 
electromagnetic wave is presented in the form: 

E = Eexp[iuJo{z/c - t)] + C.C. , 
where loq is a reference frequency and E is slowly- varying amplitude |ll|. As it was shown in 
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for a SASE PEL, driven by an electron beam with linear energy chirp, E can be written 
as follows (we use notations from |ll|): 

E = 2EoJ2 e-''^/''e^'^'^^'-^/^-'^^giz, s - %, a) (1) 

j 

Here p = Awr/(47r) is the efficiency parameter, T^ = ttJoK^Ajj/^IaXwIq), jo is the beam 
current density, Ia = mc^/e ^ 17 kA, 70 is relativistic factor, K = e\„H„ / {2\^Tcmc'^) is rms 
undulator parameter, Ajj = Jq{Q) — Ji{Q) is the Bessel function factor, Q = K^/[2{l+K'^)], 
Eq = pT'jQmc^ / {eK Ajjy/2) is the saturation field amplitude, z = F^; is a normalized position 
along the undulator, s = puJo{z/vzo — t) is normalized position along the electron bunch, 
Vzo is average longitudinal velocity (defined for a reference particle). Let the energy linearly 
depend on a particle position in the bunch (or arrival time). The energy chirp parameter 



« = -—:r—o (2) 



dt -foUJop^ 

is defined such that, for positive sign of a, particles in the head of the bunch have larger 
energy than those in the tail. Relativistic factor 70 for a reference particle (placed at s = 0) 
and reference frequency ujq are connected by the PEL resonance condition: Uq = 2c/c^7q/(1 + 
K"^). Note that the theory is applicable when pd <^ 1 Js!]. It is also useful to define normalized 



detuning 



li: C=[k^-uj{l + K^)/2c-fi]/r. 



The Green's function g, entering Eq. (JJ, is given by the inverse Laplace transform 
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exp[f{p,z,s,a)] 
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where 



f{p,z,s,a) =p{z-2s) + 
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We use a saddle point approximation to get an estimate of the integral Q for large values 
of z J8|. The saddle point is determined from the condition /' = which leads to the 4th 
power equation with three parameters: 



p'^ + 2iasp^ 



2 2 2 4zs 2ds^ 

asp— - — — p + 
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Once the saddle point, po, is found, the Green's function can be approximated as follows: 



/. . .X '2exp[f{po,z,s,a)] . 

«<^-^'"' = P„|2./>„,i,j,a)iv^ (*" 

Let us first consider the case when the energy chirp is a small perturbation, \a\z <^ 1, 
z ^ 1. A second-order expansion of the Green's function takes the following form 
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The leading correction term is the last term in the argument of the exponential function. 
This term was found in 8] (note difference in definition of normalized parameters). Setting 
d = 0, one gets from ((Zj) the well-known Green's function for unchirped beam (l2| |. 

Now let us consider the case d > and 1 ^ d ^ f . The Green's function for s ^ a^^ 
is approximated by: 



g{z,s,a) 
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More thorough analysis for small values of s shows that the Green's function has a maximum 
at Sm = 2^/^d~^, i.e. the position of maximum is independent of z while the width of 



the radiation wavepacket is proportional to A/dl. The mean frequency of the radiation 
wavepacket corresponds to a resonant frequency at s = 0. Note also that the beam density 
excitation is concentrated near s = within much shorter range, of the order of d'^/^i"^/^. 
In the case of a < and 1 ^ |a| ^ i the Green's function is given by: 

2i/4e--/2 / [^ 

g{z,S,a) ~ 1/21 -15/4-3/4- ^^P 2Wpq 
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The width of the radiation wavepacket (and of the beam density excitation as well) is of the 
order of \a\~'^^^z~^^^. The maximum of the wavepacket is positioned at Sm = 2^^^\a\~'^^^z~^^^, 
i.e. the wavepacket is shrinking and back-propagating (with respect to the electron beam) 
with increasing z. The mean frequency of the wavepacket is blue-shifted with respect to 
resonant frequency at s = 0. In normalized form this shift is AC = —\a\z/2. 

III. LINEAR REGIME OF SASE FEL 

The normalized radiation power (normalized efficiency), < fj >= -PsASE/p-Pbeam, can be 
expressed as follows |ll|: 

<.>-4^, (10) 

where < ... > means ensemble average. One can easily get from (^: 

1 f°° 

<fi{z,a) >= — ds\g{z,s,a)\'^ . (11) 

Here A^c = ^A/(2vrp) is a number of cooperating electrons (populating As = 1), A^a is a 
number of electrons per wavelength. The local power growth rate [l3| can be computed as 
follows: 

G(z,a) = --\n<fi(z,a) > . (12) 

dz 



Applying Eqs. (|TT|). (fT^ to the asymptotical cases, considered in the previous Section, 
we get the following results. For the case |d|z -C 1, z ^ 1 the FEL power is given by 
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and the local power growth rate is 



G{z,a) ~ Vs 
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It reaches maximum Gm = v^ 1 — (|d;|/16)^'^ +d/12 at the position Zm = 3^/^2^/^/|dp/^. 
Numerical simulations show that Eqs. ()13|) and ()14|) are pretty accurate up to the values 
\a\z of the order of unity although the condition \a\z <^ 1 was used to derive them. 
For the case d > and 1 <^ d <^ z we get rather simple expressions: 
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For large negative values of a we obtain: 



< fj >'- 
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IV. NONLINEAR REGIME 



We studied nonlinear regime of a chirped SASE FEL operation with 1-D version of the 



code FAST 
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Il4| . Analytical results, presented above, were used as a primary standard 
for testing the code in linear regime. Green's function was modelled by exciting density 
modulation on a short scale. As ^ 1. SASE FEL initial conditions were simulated in a 
.a„d„d wa. H. The ..1. of nu^enca, ..ulatio. . aU ca.e. we. . a good ag.ee.ent 
with analytical results presented in two previous Sections. 




FIG. 1: Increase of saturation length Azgat = -Zsatlo) — -ZsatlO) versus parameter a. Here Zsat(O) = 
13. 

The main results of the simulations of the nonlinear regime are presented in Figs. ^ and |21 
Saturation length and power are functions of two parameters, a and A'^c. For our simulations 
we have chosen A'"c = 3 x 10^ - a typical value for VUV SASE FELs. Note, however, that 
the results, presented in Figs. ^ and |H very weakly depend on A'"c. Fig. ^ shows increase of 
saturation length with respect to unchirped beam case. In Fig.|21the output power is plotted 
versus chirp parameter for two cases: when undulator length is equal to a saturation length 
for a given a and when it is equal to the saturation length for the unchirped beam case. 
One can see sharp reduction of power for negative a while a mild positive chirp [a < 0.5) is 
beneficial for SASE. 

V. ENERGY CHIRP AND UNDULATOR TAPERING 



Let us consider now the case when there is no energy chirp {a = 0) and the detuning 

I [ 
parameter changes linearly along the undulator [11]: C{z) = hiz. This change can be 

due to variation of undulator parameters {K{z) and/or A;„(z)), or due to an energy change 

7o(-S). We have found from numerical simulations that in such case the effect on FEL gain 

is exactly the same as in the case of energy chirp and no taper if d = 26i for any value 
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FIG. 2: Normalized output power versus parameter a. Solid: z = Zsat(a) (see Fig. ^; dash: 

Z = isat(O) = 13. 

of a (Figure ini shows an example). Therefore, all the results of two previous Sections can 
be also used for the case of linear variation of energy or undulator parameters with the 
substitution a — ^ 26i. The amplitudes of Green's functions are also the same while the 
phases are obviously different. In case of 6i = 0, d 7^ there is a frequency chirp along the 
bunch while in the case 61 = 0, d 7^ the frequency is changing along the undulator. 

An effect of undulator tapering (or energy change along the undulator) on FEL gain was 
studied in |lj] in the limit 61 <^ 1. Comparing our Eq. (fT!?|l (with the substitution d -^ 2h\) 
and Eq. (45) of Ref. [13|, we can see that quadratic correction term in the argument of the 



exponential function is the same but the linear term is two times larger in 



1311. The reason 



for discrepancy is that the frequency dependence of the pre-exponential factor in Eq. (42) 



screp; 
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of Ref. [13 is neglected. 

A symmetry between two considered effects (energy chirp and undulator tapering) can be 
understood as follows. If we look at the radiation field acting on some test electron from an 
electron behind it, this field was emitted at a retarded time. In the first case a back electron 
has a detuning due to an energy offset, in the second case it has the same detuning because 
undulator parameters were different at a retarded time. The question arises: can these two 




FIG. 3: Normalized power versus undulator length. Solid line 1: a = 0, bi = 0; triangles: a = 4, 
bi = —2; solid line 2: a = 4, 6i = 0; circles: a = 0, 6i = 2. 



effects compensate each other? We give a positive answer based on numerical simulations 
(see Fig. IHlas an example): by setting 6i = —a/2 we get rid of gain degradation, and FEL 
power at any point along the undulator is the same as in the case of unchirped beam and 
untapered undulator. This holds for any value of a. For instance, if one linearly changes 
magnetic field H^ of the undulator, the compensation condition can be written as follows 
(nominal values of parameters are marked with subscript '0'): 

(19) 

Of course, in such a case we get frequency chirped SASE pulse. Since compensation of gain 
degradation is possible also for large values of a (there is no theoretical limit on the value 
of chirp parameter, except for above mentioned condition pd <^ 1), one can, in principle, 
organize a regime when a frequency chirp within an intensity spike is much larger than the 
natural FEL bandwidth (given by pujo). 
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FIG. 4: Schematic diagram of attosecond X-ray source. Energy modulator performs slice energy 
modulation of the electron bunch (see Fig. ^. Undulator tapering leads to complete suppression 
of the amplification process in the most fraction of the electron bunch, and output X-ray pulse has 
200 attosecond pulse duration. 

VI. GENERATION OF ATTOSECOND PULSES 

schemes for gen eration of attosecond pulses from X-ray SASE FELs 
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Up to now several 
have been proposed 

llSJ I making use of energy modulation of a short slice in the electron bunch by a high-power 
few-cycle optical pulse in a two-period undulator. Due to energy modulation the frequency 
of SASE radiation in X-ray undulator is correlated to the longitudinal position within the 
few-cycle-driven slice of the electron beam. The largest frequency offset corresponds to a 
single-spike pulse in time domain (about 300 attoseconds). The selection of single-spike 
pulses is achieved by using a crystal monochromator after the X-ray undulator [l7|, or with 
the help of the other undulator tuned to the offset frequency [la| • 

In this paper we propose a new scheme (see Fig. |3]) that makes use of the compensation 
effect, described in the previous Section. Indeed, there is a strong energy chirp around zero- 
crossing of energy modulation. If one uses appropriate undulator taper then only a short 
slice around zero-crossing produces powerful FEL pulse. The main part of the bunch is 
unmodulated and suffers from strong negative undulator tapering (see Fig. |21). One should 
also note that for large negative taper the SASE FEL gain is very sensitive to longitudinal 
velocity spread. Therefore, a high-contrast attosecond pulse is directly produced in the 
undulator. 
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FIG. 5: Possible evolutions of the electric field in the 5-fs pulse, carried at a wavelength 800 nm 
for two different pulse phases {4> = 0, 7r/2) 



Operation of attosecond SASE FEL is illustrated for the parameters close to those of the 
European XFEL operating at the wavelength 0.15 nm 9|. The parameters of the electron 
beam are: energy 15 GeV, charge 1 nC, rms pulse length 25 fim, rms normalized emittance 
1.4 mm-mrad, rms energy spread 1 MeV. Undulator period is 3.65 cm. 

The parameters of the seed laser are: wavelength 800 nm, energy in the laser pulse 3 mJ, 
and FWHM pulse duration 5 fs (see Fig. \^. The laser beam is focused onto the electron 
beam in a short undulator resonant at the optical wavelength of 800 nm. Optimal conditions 
of the focusing correspond to the positioning of the laser beam waist in the center of the 
modulator undulator. It is assumed that the phase of laser field corresponds to "sine" mode 
(dashed line with (p = tt/2, see Fig. El). Parameters of the modulator undulator are: period 
length 50 cm, peak field 1.6 T, number of periods 2. The interaction with the laser light 
in the undulator produces a time-dependent electron energy modulation as it is shown in 
Fig. ini This modulation corresponds to the energy chirp parameter d ^ 2 at zero crossing 
(t = 5fsinFig.ini). 

Optimization of the attosecond SASE FEL has been performed with the three- 
dimensional, time dependent code FAST jl4| taking into account all physical effects in- 
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FIG. 6: Energy modulation of the electron beam at the exit of the modulator undulator (dotted 
line) and a profile of the radiation pulse at the undulator length 100 m 

fluencing the SASE FEL operation (diffraction effects, energy spread, emittance, slippage 
effect, etc.). Three-dimensional simulations confirmed the predictions of the one-dimensional 
model: the energy chirp and the undulator tapering compensate each other, there is strong 
suppression of the amplification in the case of uncompensated negative taper. 

Undulator tapering is performed by changing the gap of undulator modules [2| such that 
magnetic field increases linearly along the undulator length (6i < 0). We performed the 
scan of tapering depth bi in order to maximize the power in the main peak on one hand, 
and to minimize contribution of the background, on the other hand. We ended up with the 
value of taper which is about 20 % smaller than that required for a perfect compensation of 
chirp at t = 5 fs. Note that the chirp is not linear in the region of interest. In addition, a 
mild net positive chirp is beneficial for SASE, as it was discussed above (see Fig. |2I). 

A typical radiation pulse at the undulator length 100 m is shown in Fig. IHl One can see 
a high-power spike in the region where the energy chirp is well compensated by the taper 
and two weak side peaks at t ~ 2 fs and t ~ 8 f s where the net effect is negative taper. In 
the rest of the bunch a large negative taper together with velocity spread and 3-D effects 
completely suppresses amplification. In Fig. [7| we present three different shots illustrating 
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FIG. 7: Temporal structure of the radiation pulse (three different shots) at the undulator length 
100 m. 

the properties of the main peak. Typical pulse duration is about 200 attoseconds (FWHM) 
and peak power ranges from several tens up to hundred GW. To estimate the contrast 
(which we define as the ratio of energy in the main peak to the total radiated energy at the 
experiment) we assume that an angular coUimation is used in order to reduce spontaneous 
emission background. A collimator with half-angle 3 /irad allows the entire intensity in the 
main peak to be transmitted. The contrast is influenced by SASE intensity in two side 
peaks and by spontaneous emission in the first harmonic from the rest of the bunch. For the 
charge of 1 nC, as in our numerical example, the contrast is about 95 %. Higher harmonics of 
undulator radiation (if they disturb an experiment) can be cut, for instance, by a multilayer 
monochromator with a bandwidth of the order of 1 %. 

VII. BEYOND "FUNDAMENTAL LIMIT" 



It is generally accepted that the shortest pulse, that can be obtained from a SASE FEL, 
is given by a duration of intensity spike in time domain, i.e. it is defined by inverse FEL 
bandwidth (pcjo)^^- However, the fact that a SASE FEL can operate with a strong chirp 
parameter (in combination with undulator tapering) without gain degradation, opens up a 
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possibility of a conceptual breakthrough: one can get from SASE FEL a radiation pulse 
which is much shorter than the inverse FEL bandwidth. Indeed, in the case of d ^ 1, 
the frequency chirp inside an intensity spike is much larger than FEL bandwidth. Thus, 
one can use a monochromator to reduce pulse duration. By an appropriate choice of the 
monochromator bandwidth one can select an X-ray pulse that is shorter by a factor of 
v2A than the inverse FEL bandwidth. The only theoretical limit in this case is given by 
the condition pa <^ 1. Note that for hard X-ray FELs the parameter p is in the range 
10-^ - 10-3. 

To illustrate a possible technical realization of this idea, we can suppose that the energy 
modulation by a few-cycle optical pulse is increased by a factor 3 so that d ~ 6. In 
combination with undulator tapering and a monochromator, this would allow to obtain sub- 
100-GW coherent X-ray pulses with a duration below 100 attoseconds and a contrast about 
80-90 %. 
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